Control of unstable steady states by extended time-delayed feedback 



O 

o 

< 



U 
d 



> 



oo 

o 

O 



X 



Thomas Dahms, Philipp Hovel, and Eckehard ScholEI 

Institut fur Theoretische Physik, Technische Universitdt Berlin, 10623 Berlin, Germany 

(Dated: February 1, 2008) 

Time-delayed feedback methods can be used to control unstable periodic orbits as well as unstable 
steady states. We present an application of extended time delay autosynchronization introduced by 
Socolar et al. to an unstable focus. This system represents a generic model of an unstable steady 
state which can be found for instance in a Hopf bifurcation. In addition to the original controller 
design, we investigate effects of control loop latency and a bandpass filter on the domain of control. 
Furthermore, we consider coupling of the control force to the system via a rotational coupling matrix 
parametrized by a variable phase. We present an analysis of the domain of control and support our 
results by numerical calculations. 
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I. INTRODUCTION 

Since the last decade, the stabilization of unstable and 
chaotic systems has been a field of extensive research. A 
variety of control schemes have been developed to control 
periodic orbits as well as steady states [H, [2, S. A simple 
and efficient scheme, introduced by Pyragas [4], is known 
as time delay autosynchronization (TDAS). This control 
method generates a feedback from the difference of the 
current state of a system to its counterpart some time 
units r in the past. Thus, the control scheme does not 
rely on a reference system and has only a small number 
of control parameters, i.e., the feedback gain K and time 
delay r. It has been shown that TDAS can stabilize both 
unstableperiodic orbits, e.g., embedded in a strange at- 
tractor jj [Bj] , and unstable steady states [1, 0, • In the 
first case, TDAS is most efficient if r corresponds to an 
integer multiple of the minimal period of the orbit. In 
the latter case, the method works best if the time delay 
is related to an intrinsic characteristic timescale given by 
the imaginary part of the system's eigenvalue Q. A gen- 
eralization of the original Pyragas scheme, suggested by 
Socolar et al. [9], uses multiple time delays. This ex- 
tended time delay autosynchronization (ETDAS) intro- 
duces a memory parameter R, which serves as a weight 
of states further in the past. A variety of analytic re- 
sults about time-delayed feedback control are also known 
EH 03 ' f° r instance, in the case of long time de- 
or the odd number limitation [15[, which was 
refuted recently [lr| . 

Although there has been strong effort on the research 
on the original Pyragas method [l?], EH, much less is 
known in the case of extended time-delayed feedback 
EE HI EI SBQ. Recently it was shown that 
the additional memory parameter introduces a second 
timescale which leads to a strong improvement of the sta- 
bilization ability, for instance, arbitrary large correlations 




of stochastic oscillations without inducing a bifurcation 

In the present paper, we apply the ETDAS control 
method, which was initially invented to stabilize periodic 
orbits, to an unstable steady state realized as an unsta- 
ble focus. This can be seen as a system close to, but 
above, a supercritical Hopf bifurcation. As a modifica- 
tion, we consider also an additional control loop latency, 
a bandpass filter, and different couplings. 

This paper is organized as follows: In Sec. [Ill we intro- 
duce our model equations and develop the analytic tools 
used throughout the paper. Section HIJ deals with a non- 
diagonal coupling implemented by a rotational matrix. 
In Sec. IIVI we consider latency time effects, which arise 
if a time lag exists between calculation of the control force 
and reinjection into the system. Section [V] introduces a 
specific modification of ETDAS that includes a bandpass 
filtering of the control signal. This is important if high 
frequency components are present in the system. Finally, 
we conclude with Sec. EH 



II. EXTENDED TIME-DELAYED FEEDBACK 

We consider an unstable fixed point of focus type. 
Without loss of generality, the fixed point is located at 
the origin. In complex center manifold coordinates z the 
linearized system can be written as z(t) = (A + iw)z(i), 
where A and uo are real numbers corresponding to the 
damping and oscillation frequency, respectively. The sta- 
bility of the steady state is determined by the sign of the 
real part of the complex eigenvalue A + iuo. Since we 
consider an unstable focus, e.g., closely above a Hopf bi- 
furcation, we choose the parameter A to be positive and 
uj nonzero. Separated in real and imaginary parts, i.e., 
z(t) = x(t) + iy{t), the dynamics of the system is given 

by 



±{t) = A x(i) - F(t), 



(1) 
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where x is the state vector composed of the real and 
imaginary part x and y of the variable z, the matrix A 
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denotes the dynamics of the uncontrolled system, 

a = (-Ia). < 2 > 

and F is the ETDAS control force, which can be written 
in three equivalent forms 

oo 

F(t) = K Rn I x (* " nT ) ~ x (* " ( n + i) 7 ")] ( 3 ) 



K 



x(t) - (1 - i?) ^ iT"^ - nr) 



if [x(t) - x(t - t)] + iJF(t - r), 



(4) 
(5) 



where K and r denote the (real) feedback gain and the 
time delay, respectively. R G (—1, 1) is a memory param- 
eter that takes into account those states that are delayed 
by more than one time interval r. Note that R = yields 
the TDAS control scheme introduced by Pyragas 

The control force applied to the ith component of the 
system consists only of contributions of the same compo- 
nent. Thus, this realization is called diagonal coupling. 
We will consider a generalization to a nondiagonal cou- 
pling scheme in Sec. IIIII The first form of the control 
force, Eq. (j3]), indicates the noninvasiveness of the ET- 
DAS method because x*(t — r) = x*(t) if the fixed point 
is stabilized. The third form, Eq. ([5]), is suited best for 
an experimental implementation since it involves states 
further than r in the past only recursively. 

While the stability of the fixed point in the absence of 
control is given by the eigenvalues of matrix A, i.e., A ± 
iuj, one has to solve the following characteristic equation 
in the case of an ETDAS control force: 



-At 



k + K 



1 - Re~ Ar 



(6) 



Due to the presence of the time delay r, this character- 
istic equation becomes transcendental and possesses an 
infinite but countable set of complex solutions A. In the 
case of TDAS, i.e., R = 0, the characteristic equation can 
be solved analytically in terms of the Lambert function 
d, [H, [27], [2^ • We stress that for nonzero memory pa- 
rameter R, however, such a compact analytic expression 
is not possible. Thus, one has to solve Eq. (j6j) numeri- 
cally. 

Figure [I] depicts the dependence of the largest real 
parts of the eigenvalue A upon the time delay r ac- 
cording to Eq. (|6]) for different memory parameters R 
and fixed feedback gain K = 0.3. The red, green, 
black, blue, and magenta curves of Re(A) correspond to 
R = —0.7,-0.35,0,0.35,0.7, respectively. The param- 
eters of the unstable focus are chosen as A = 0.1 and 
ou = 7r. Note that the time delay r is given in units 
of the intrinsic period To = 2tt/oj. When no control is 
applied to the system, i.e., r = 0, all curves start at A 
which corresponds to the real part of the uncontrolled 
eigenvalue. For increasing time delay, the real part of A 




FIG. 1: (Color online) Largest real part of the complex eigen- 
values A as a function of r for different values of R. The 
red, green, black, blue, and magenta curves correspond to 
R = —0.7, —0.35, 0, 0.35, 0.7, respectively. The parameters of 
the unstable focus are chosen as A = 0.1 and uj — tt which 
yields an intrinsic period To = 27t/lj = 2. The feedback gain 
K is fixed at KT = 0.6. 



decreases and eventually changes sign. Thus, the fixed 
point becomes stable. Note that there is a minimum of 
Re (A) indicating strongest stability if the time delay r 
is equal to half the intrinsic period. For larger values 
of r, the real part increases and becomes positive again. 
Hence, the system loses its stability. Above r = To, the 
cycle is repeated but the minimum of Re (A) is not so 
deep. The control method is less effective because the 
system has already evolved further away from the fixed 
point. For vanishing memory parameter R = (TDAS), 
the minimum is deepest, however, the control interval, 
i.e., values of r with negative real parts of A, increases 
for larger R. Therefore the ETDAS control method is 
superior in comparison to the Pyragas scheme. 

Figure [2] shows the domain of control in the plane 
parametrized by the feedback gain K and time delay 
r for different values of R : 0,0.35,0.7,-0.35 in pan- 
els (a), (b), (c), and (d), respectively. The color code 
indicates only negative values of the largest real parts 
of the complex eigenvalue A. Therefore, Fig. [1] can be 
understood as a vertical cut through Fig. [2] for a fixed 
value of KTo = 0.6. Each panel displays several islands 
of stability which shrink for larger time delays r. Note 
that no stabilization is possible if r is equal to an integer 
multiple of the intrinsic period To. The domains of con- 
trol become larger if the memory parameter R is closer 
to 1. 

In order to obtain some analytic information of the do- 
main of control, it is helpful to separate the characteristic 
equation (j6]) into real and imaginary parts. This yields 
using A = p + iq: 



K(l -e~ pr cos qr) = A - p - Re~ 



:[(A — p) cos qr ± (uj — q) sin qr] (7) 



3 



1 1 1 1 , J, 1 
(a) R = 


1 1 1 1 „ J r, 1 1 
^^^^^^^^^^ __ IUI 1 1 — U.JJ 










_ i i i i i i i 
_BMm^^^^^R_= 0.7_ 


1 I I I 1 1 1 

(d) K = -0.35_ 













-0.2 
-0.4 
-0.6 
-0.8 



-0.2 
-0.4 
-0.6 
-0.8 
-1 



2 4 6 
KT n 



2 4 

KT n 



FIG. 2: (Color online) Domain of control in the (K,r) plane 
for different values of R : 0, 0.35, 0.7, —0.35 in panels (a), (b), 
(c), and (d), respectively. The greyscale (color online) shows 
only negative values of the largest real part of the complex 
eigenvalues A according to Eq. JB]). The parameters of the 
system are as in Fig. [T] 



and 



Ke~ pr sin qr = ±(<j - q) + Re~ pr 

x [(A — p) sin qr db (co — q) cos qr] . (8) 

The boundary of the domain of controls is determined 
by a vanishing real part of A, i.e., p = 0. With this 
constraint, Eqs. ([7j) and (j8|) can be rewritten as 

K(l — cos qr) = A — R[\ cos qr ± (a; — g) sin gr], (9) 
ifsingr = ±(a; — q) + ir![Asingr ± (a; — q) cos gr]. 

At the threshold of control (p = 0,<? = a;), there is 
a certain value of the time delay, which will serve as a 
reference in the following Section, given by 



(2n + 1)tt 



To 



(10) 



where n is any nonnegative integer. For this special 
choice of the time delay, the range of possible feedback 
gains K in the domain of control becomes largest as can 
be seen in Fig. [2j Hence, we will refer to this r value 
as optimal time delay in the following. The minimum 
feedback gain at this r can be obtained: 



-K-min(R) — 



\(1 + R) 



(11) 



Extracting an expression for sin(gr) from Eq. ([9|) and 
inserting it into the equation for the imaginary part leads 
after some algebraic manipulation to a general depen- 
dence of K on the imaginary part q of A 



K(q) = 



(1 + R) [A 2 + (a, - g) 2 ] 



2A 



(12) 



Taking into account the multivalued properties of the 
arcsine function, this yields in turn analytical expressions 
of the time delay in dependence on q 



ri(q) = 



T 2{q) 



arcsm 



2A(l-ir)(a;-g) 

A 2(l_ jR 2)2 + ( u; _ (? )2( 1 + jR )2 



2nir 



Q 

• ( 2X(l-R 2 )(uj-q) 

•arcsm ^ A 2( 1 _ jR2 ) 2+ ( a; _ g )2( 1+jR )2 



(2n- 



(13) 
1)tt 



where n is a nonnegative integer. Together with Eq. ([T2]h 
these formulas describe the boundary of the domain of 
control in Fig. [2j Note that two expressions r\ and T2 are 
necessary to capture the complete boundary. The case 
of TDAS control was analyzed in [14] and is included as 
special choice of R = 0. 




FIG. 3: (Color online) Domain of control in the (K, R) plane 
for different values of A. The blue, red, green, and yellow 
domains correspond to ATb = 0.2,1,5, and 10, respectively, 
as indicated. The time delay is chosen as r = To/2 and uo — tt. 

For a better understanding of effects due to the mem- 
ory parameter R, it is instructive to consider the domain 
of control in the plane parametrized by R and the feed- 
back gain K. The results can be seen in Fig. [3l where 
the blue, red, green, and yellow areas correspond to the 
domain of control for ATb = 0.2, 1, 5, and 10, respectively. 
The other system parameter is chosen as uo = tt. We keep 
the time delay constant at r = To/2. Note that the K 
interval for successful control increases for larger values 
of R. In fact, while the original Pyragas scheme, i.e., 
R = 0, fails for AT = 10, the ETDAS method is still 
able to stabilize the fixed point. The upper left bound- 
ary corresponds to Eq. ([TT]) . The lower right boundary 
can be described by a parametric representation which 
can be derived from the characteristic equation ([6]): 



R 



Kr 



Ar - tftan (tf/2) 
Ar + tftan (tf/2) ' 

ti 2 + (Ar) 2 
Ar + tftan (tf/2) ' 



(14) 
(15) 



where we used the abbreviation $ = (q — uj) r for nota- 
tional convenience. The range of d is given by d G [0, tt). 
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A linear approximation leads to an analytic dependence 
of R an the feedback gain K given by a function R(K) 
instead of the parametric equations ([T4]) and ([T5|) . A 
Taylor expansion around $ = tt yields 



-K-max (R~) 



A 2 +- 



2A 



(R+1) + 2(R-1). 



(16) 




K T n 



FIG. 4: (Color online) Domain of control in the (K, A) plane 
for different memory parameters R. The yellow, red, green, 
and blue areas correspond to R — —0.35, (TDAS), 0.35, 
and 0.7, respectively. The time delay is fixed at r = To/2. 

Another representation of the superior control ability 
of ETDAS is depicted in Fig. [4j The domain of control 
is given in the {K, A) plan for different values of R. The 
yellow, red, green, and blue areas refer to R = —0.35, 
(TDAS), 0.35, and 0.7, respectively. The time delay is 
chosen as r = To/2. One can see that for increasing i?, 
the ETDAS method can stabilize systems in a larger A 
range. However, the corresponding K interval for suc- 
cessful control can become small. See, for instance, the 
blue area (R = 0.7) for large A. A similar behavior was 
found in the case of stabilization of an unstable periodic 
orbit by ETDAS j29j]. We stress that, as in the case of 
periodic orbits, the boundaries of the shaded areas can 
be calculated analytically from the following expression: 



Kt 



At 



(1 - Ryd f l + R 

tan(tf/2) V 1 -^ 

t? f l + R 

tan(tf/2) \1-R 



-tan 2 (tf/2) 



(17) 



(18) 



where we used $ = (q — uo)r with d G [0, tt) as in 
Eqs. (p3j) and (fT5|). The maximum value for A, which 
can be stabilized, is given by the special case $ = 0: 



Xn 



cT = 2 



R 



1-R 



(19) 



In this Section, we have focussed our discussion to the 
ETDAS method in the simplest realization of diagonal 
coupling. In the following, we will investigate the effects 
of nondiagonal coupling introduced by a variable phase. 



III. PHASE DEPENDENT COUPLING 

Time-delayed feedback has been widely used in opti- 
cal systems both to study the intrinsic dynamics and to 
control the stabilization of, for instance, a laser device 
[30l l3ll I32I l33l [H, [35[. In these systems, the optical 
phase is an additional degree of freedom. We consider 
this additional control parameter as a generalization of 
the ETDAS feedback scheme using a nondiagonal cou- 
pling as opposed to the diagonal coupling discussed in the 
previous Section. This nondiagonal coupling is realized 
by introducing a coupling matrix containing a variable 
phase up: 



A uo 
—uo A 

cos up 
sm up 



x 

y 

— sin up 
cos up 



F(t), 



(20) 



In optical systems like semiconductor lasers with ex- 
ternal optical feedback [36|, [33] , this feedback phase can 
be seen as the phase of the electric field. Experimen- 
tally, this phase of the feedback can be varied by tuning 
an external Fabry-Perot cavity. It has also been demon- 
strated that a feedback phase plays an important role 
in the suppression of collective synchrony in a globally 
coupled oscillator network |38| . 




FIG. 5: (Color online) Largest real part of the eigenvalues 
A as a function of r for different phases ip. The black, red, 
green, and blue curves correspond to ip = 0, 7r/8, tt/4, and 
37r/8, respectively. The other control parameters are fixed as 
R = 0.7 and KTo =0.6. The parameters of the system are 
as in Fig. [1] 

The stability of the fixed point is again given by the 
largest real part of the complex eigenvalues A, which are 
calculated as the solutions of the following modified char- 
acteristic equation: 



A + Ke^ 



1 



-At 



1 - Re~ Ar 



A + iuo. 



(21) 
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FIG. 6: (Color online) Domain of control in the (cp, il~) plane 
for optimal time delay r = To/2. Panels (a), (b), (c), and 
(d) correspond to a memory parameter R of 0, 0.35, 0.7, and 
—0.35, respectively. The color code shows the largest real 
part of the complex eigenvalues A as given by Eq. (|2T|) . The 
parameters of the system are as in Fig. [1] 
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(22) 



Note that the time delay that corresponds to this value 
of K m i n is no longer given by Eq. (jTUJ) and is not the 
optimal time delay in the general case of nonzero phase. 
Nevertheless, Eq. ([22]) can be used as a coarse estimate 
of the minimum feedback gain for the regime of small 
values of cp, if the time delay is chosen as Eq. (fit 



-k/2 Jt/2 
cp 
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FIG. 7: (Color online) Domain of control in the (cp, K) plane 
for time delay r = O.lTo. Panels (a), (b), (c), and (d) corre- 
spond to a memory parameter R of 0, 0.35, 0.7, and —0.35, 
respectively. The color code shows the largest real part of the 
complex eigenvalues A as given by Eq. (|2T|) . Only negative 
values are depicted. The parameters of the system are as in 

Fig.m 



Note that this equation differs from the characteristic 
equation in the diagonal case (see Sec. [TTj) by an addi- 
tional exponential term. This is due to the choice of 
phase dependent coupling by a rotational matrix. We 
stress that a similar phase factor has recently been used 
[lq to overcome some topological limitation of time- 
delayed feedback control known as odd number limita- 
tion theorem, which refers to the case of an unstable peri- 
odic orbit with an odd number of real Floquet multipliers 
larger than unity. [13, [39[ 

In analogy to Sec. [Til a minimum feedback gain can be 



FIG. 8: (Color online) Domain of control in the (K, R) plane 
for different values of the feedback phase ip. The blue, green, 
red, and yellow areas correspond to ip = 0, 7r/8, tt/4, and 
37r/8, respectively. Panel (a) displays the domain of control 
for optimal r = To/2; panel (b) for r = 0.1Tb. The parame- 
ters of the system are as in Fig. [T] 

Figure [5] depicts the dependence of the largest real part 
of the eigenvalues A on the time delay r for fixed values 
of R = 0.7 and KTq = 0.6, but different values of the 
phase. The black, red, green, and blue curves correspond 
to ip = 0, 7r/8, 7r/4, and 37r/8, respectively. It can be 
observed that the control is overall less effective for larger 
cp, as the curves are shifted up towards positive reals parts 
for increasing the phase. The range of possible values for 
the time delay shrinks. The optimal time delay is shifted 
towards smaller values for larger cp, which can be seen for 
the case of cp = 37r/8, where the optimal time delay is in 
the range of r = O.ITq instead of 0.5Tq, which was the 
optimal time delay for cp = according to Eq. ([T0|) . 
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The four-dimensional parameter space is now given by 
the feedback gain time delay r, memory parameter 
and feedback phase tp. At first, we consider the domain 
of control in the plane parametrized by K and cp. Hence, 
we keep the other remaining control parameters R and cp 
fixed. Figures [6] and show the domain of control for a 
time delay of To/2 and 0.1Tb, respectively. In each figure, 
the memory parameter R is chosen as R = 0, 0.35, 0.7, 
and —0.35 in panels (a), (b), (c), and (d), respectively. 
The color code corresponds to the largest real part of the 
complex eigenvalues as calculated from Eq. ([2T]) . Only 
negative values are depicted, i.e., those combinations of 
K and cp for which the control scheme is successful. Note 
that the case R = corresponds to the TDAS control 
method [37} • An increase of the memory parameter R 
leads to a larger domain of control. Even though the 
system can be stabilized for a larger range of K and cp, 
the system becomes over all less stable, since the real 
part of A is closer to zero. For negative values of R, the 



domain of control shrinks. Note that also in the case of 
non-optimal time delay as in Fig. [71 the range of choice 
for possible feedback gain and phase is enlarged. 

For a better understanding of the effects of the feed- 
back phase, Fig. [8] depicts the domain of control in the 
(K, R) plane. The blue, green, red, and yellow areas cor- 
respond to successful control for cp = 0, 7r/8, 7r/4, and 
37r/8, respectively. Panel (a) shows the case of optimal 
time delay, i.e., r = To/2; panel (b) displays the case of 
r = O.lTo. Note that an increase of cp leads to a smaller 
domain of control in the case of r = To/2. This effect, 
however, is reversed for non-optimal choices of r, where 
the phase cp compensates for the bad choice of the time 
delay. Thus, control is possible again, for instance, in the 
TDAS case (R = 0). Following the strategy introduced 
in Sec. HI one can derive also in the case cp ^ para- 
metric formulas for the the boundary of the domain of 
control: 



R 
Kt 



d [cos (cor + $ + op) — cos op] + At [sin op — sin (out + $ + op)] 
d [cos cp — cos (our + d + — At [sin cp + sin (our + $ + (/?)]' 

(tf 2 + A 2 T 2 )cOS [KwT + 0)] 

At cos [\ (cor + #) - cp] - sin [\ (cor + 



(23) 
(24) 



We stress that it is possible to derive expression of 
K(q) and r(q) similar to Eqs. ([T2]) and (jTTJ) also in the 
case of (p 7^ 0. These calculations are lengthy and do not 
produce more insight and thus are omitted here. 



IV. CONTROL LOOP LATENCY 

After the investigation of nondiagonal coupling, we 
consider in this Section an additional control loop la- 
tency. This latency is associated with time required for 
the generation of the control force and its reinject ion into 
the system. In an optical realization, for instance, the la- 
tency time is given by the propagation time of the light 
between the laser and the Fabry-Perot control device. We 
stress that in the case of unstable periodic orbits, Just 
has shown that longer latency times reduce the control 
abilities of the time-delayed feedback of TDAS type [29[ . 
Similar results were found for ETDAS [H[ . 

The latency time 5 acts as an additional time delay in 
all arguments of the control force of Eq. ([3]). Using the 
recursive form of F as given by Eq. ([5|) , this yields 

F(t) = [x(t-5)-x(t-5-T)]+KF(t-T). (25) 

The characteristic equation (|6|) is now modified by an 



additional exponential factor 

A±i„ = A + Ke-** *-^ . (26) 

In contrast to the previous Section, this exponential term 
depends on the eigenvalue A itself. 

Figure [9] depicts the dependence of the largest real part 
of the eigenvalues A on the time delay r for fixed values of 
R = 0.7 and KTq = 0.6, but different latency times. The 
black, red, green, and blue curves correspond to S = 0, 
0.1, 0.3, and 0.5, respectively. It can be seen that the 
control scheme is less successful for longer latency times. 
The r interval with negative real parts of A becomes 
smaller. In the case of S = 0.25To, for instance, control 
can only be achieved in a narrow range of small r and the 
second minimum does not reach down to negative Re (A) 
anymore. In addition, the minima of the real parts are 
distorted and shifted towards smaller time delays. 

Taking also a varying feedback gain K into account, 
the domain of control can be seen in FigfTOl The other 
control parameters are fixed at r = To/2 and R — 0.7. 
Panels (a), (b), (c), and (d) correspond to values of S = 0, 
O.lr, 0.2t, and 0.3r, respectively. As in Figs. [6] and [3 of 
the previous Section, the color code corresponds to the 
largest real part of the complex eigenvalues which are 
calculated from Eq. ([25]) . Note that only negative values 
are depicted. For increasing latency time, the domains 




which is consistent with the TDAS case investigated in 



FIG. 9: (Color online) Largest real part of the eigenvalues A 
as a function of r for different latency times 8. The black, red, 
green, and blue curves correspond to S = 0, 0.05Tb, 0.15Xb, 
and 0.25Tb, respectively. The other control parameters are 
fixed as R = 0.7 and KTo = 0.6. The parameters of the 
system are as in Fig. [1] 
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FIG. 11: (Color online) Domain of control in the (K, R) plane 
for different values of the latency time S. The blue, green, red, 
and yellow areas refer to values of 6 = 0, 0.05To, O.lTo, and 
0.15To, respectively. Panel (a) corresponds to an optimal time 
delay r = To/2, panel (b) to r = To/8. The parameters of 
the system are as in Fig. [1] 




FIG. 10: (Color online) Domain of control in the (r, K) plane 
for different values of the latency time S and fixed R = 0.7. 
Panels (a), (b), (c), and (d) correspond to values of S = 0, 
0.05To, O.lTo, and 0.15To, respectively, where the time delay 
is fixed at r = To/2. The parameters of the system are as in 

Fig.m 

of control shrink. Similar to the discussion in the one- 
dimensional projection of Fig. [9] the islands are distorted 
towards smaller time delays. 

Separating the characteristic equation ([26]) into real 
and imaginary part, one can derive, in analogy to Sec. HI 
an expression for the minimum feedback gain 

Kmin(6) = J^TTVT, v W 
cos[(2n + 1)ho/t\ 



As another two-dimensional projection of the parame- 
ter space, Fig. [11] displays the domain of control in the 
(K,R) plane for different values of the latency time S. 
The blue, green, red, and yellow areas refer to values of 
6 = 0, 0.05T , O.ITq, and 0.15T , respectively. Similar to 
Fig. [8] panel (a) shows the case of optimal choice of the 
time delay r = To/2 and panel (b) refers to r = To/8. In 
the first case, the domain of control shrinks considerably 
for increasing 5, whereas in the latter case, this change 
is less pronounced. 

Similar to the previous Sections, it is possible to derive 
a parametric expression for the boundary of the domain 
of control [R(<3) and K (<$)}: 
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R = 



Kt = 



-0 cos [(lot + 0) f ] + cos [(cor + 0) (l + f )] + Ar (sin [(cjt + (9) f ] - sin [(out + (9) (l + f )] ) 
6>cos [(wt + 0) f] - (9 cos [(wt + (9) (l - I )] - Ar (sin [(ujt + 6) f ] + sin [(wr + J) (l - f )]) ' 

' + (Ar) 2 l cos [| (wr + 9)) 



Ar cos [(I - I) (wr + 0)] + sin [(£ - §) (c^r + 0)] " 



(28) 



(29) 



V. BANDPASS FILTERING 

In the context of semiconductor lasers, it has been 
shown both theoretically and experimentally that de- 
layed, bandpass filtered optical feedback is able to sur 



press the undamping of relaxation oscillations 36. [4j 
Recently, filtering of an optical feedback signal [4l| was 
also investigated for laser models of Lang-Kobayashi type 
[42| . In addition, filtered feedback has also been proven 
important in the investigation of a Hopf bifurcation [43[ , 
and in the stabilization of unstable periodic orbits by 
ETDAS in semiconductor superlattices [4J]. In order to 
model this type of modification of the control force, a 
bandpass filter acting on the feedback force is introduced 
as a Lorentzian in the frequency domain with the transfer 
function 



1 



1 



(30) 



where c^o denotes the peak of the transfer function and 
7 is half the full width at half maximum of the function. 
To introduce the filter into the system in the time do- 
main, one can add two additional differential equations 
to Eq.(pQ) such that the original two-dimensional system 
becomes four-dimensional: 



x(t) = Ax(t) - F(x(t)), 

i(t) = j(x(t)-x(t))-w y(t), 
i{t) = i{y{t)-y(t))+co x(t), 



(31) 



where x and y denote the filtered versions of x and 
respectively. Note that the feedback force F(x(t)) is now 
generated from the filtered state vector x(t) which con- 
sists of the two filtered variables x and y: 

oo 

F(x(t)) = K^R n [x(t - nr) - x(t - (n + l)r)] . (32) 

n=0 

Equivalently to the additional differential equations in 
Eq. ([3T]h the filtering, for instance in the case of a low 
pass filter, z.e., ojq = 0, can be expressed by convolution 
integrals, where the filtered counterparts of x(t) and y(t) 
are given as follows: 



x(t) 

m 
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x^e-^-^dt' (33) 
y^'y-^-^dt'. (34) 



The system of differential equations (|3Tj) yields a char- 
acteristic equation of the form 



±i [uj (A - A) + uj (7 + A)] 
1 - e~ Ar 



(35) 



1 - Re~ Ar 



+ u u - (A - A) (7 + A) . 



The minimum feedback gain can be calculated in de- 
pendence on ojq and 7 in similarity to Eq. (fTTj): 



Kmin (^0,7) 



\(1 + R) 



(c^o + uj) 



2ujquj\ 

7 



(7-Ar 



(36) 

Note that at this special value of the feedback gain K, 
the time delay is no longer given by the intrinsic value 
r = tt/uj, neither is it the optimal time delay in the gen- 
eral case of nonzero ojq and finite values of 7. Adjusting 
the time delay according to the choice of c^o and 7 leads 
to boundaries of the domain of control that can only be 
understood in the scope of a parametric representation 
of the boundaries in the (K, i?)-plane. Such parametric 
equations can be derived in analogy to Eq. ([T5]). How- 
ever, the resulting equations for the case of the filtered 
system are not shown here due to the complexity of the 
terms. 
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FIG. 12: (Color online) Largest real part of the eigenvalues 
A as a function of r for different filter widths 7 and fixed 
coo = 0. The black, red, green, and blue curves correspond 
to 7T0 = 2000, 40, 10, and 2, respectively. The other con- 
trol parameters are fixed as R = 0.7 and KTo = 0.6. The 
parameters of the system are as in Fig. [1] 
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FIG. 13: (Color online) Largest real part of the eigenvalues 
A as a function of r for different mean values of the filter ujo 
and fixed 7T0 = 10. The black, red, green, and blue curves 
correspond to ljqTq = 0, 7r, 27t, and 47T, respectively. The 
other control parameters are fixed as R = 0.7 and KTo =0.6. 
The parameters of the system are as in Fig. [1] 



Figures [12] and [13] depicts numerical solutions of the 
charactistic equation Eq. ([35]) for different choices of the 
parameters 7 and uoq. In Fig. [12] the largest real part of 
the eigenvalues is shown in dependence on the time delay 
for different values of the filter width 7 in dependence on 
the time delay. The parameter uoq is fixed to 0, which 
correspnds to the case of a low pass filter. The black, 
red, green, and blue curves correspond to 7T0 = 2000, 
40, 10, and 2, respectively. The other control parameters 
are fixed as R = 0.7 and KTo = 0.6. It can be seen, 
that large values of 7 show similar behavior as in the 
unfiltered system in Fig. [TJ Decreasing 7 flattens the 
curves and shifts them up towards positive real parts. 
Therefore, no control is possible for the case of 7X0 = 2. 
Additionally, the optimal time delay is shifted to smaller 
values, which can be seen from the case 7X0 = 10. 

In Fig. [13] similar curves are shown for fixed 7X0 = 10 
and different choices of the filter's mean value o;o. The 
black, red, green, and blue curves correspond to c^oTo = 
0, 7r, 27r, and 47r, respectively. The other control param- 
eters are fixed as R = 0.7 and KTq = 0.6. It can be 
observed that the increase of uoq shifts the curves further 
upwards to positive real parts, leading to a less stable sys- 
tem. The range of the time delay for successful control 
shrinks. 

The domain of control can also be investigated in the 
(if, R) plane. Results are shown in Fig. [TH for different 
values of the filter width 7 and fixed uoq = 0. The blue, 
green, red, and yellow areas refer to values of 7X0 = 2000, 
40, 10, and 2, respectively. The other control parameters 
are fixed as R = 0.7 and r = To/2. The domain of control 
for large values of 7, here for instance 7T0 = 2000, looks 
very similar to that in the case of the unfiltered system, 
that is depicted in Fig. [3] Decreasing the value of 7, the 
region of control shrinks. The lower right boundary is 




FIG. 14: (Color online) Domain of control in the (if, R) plane 
for different values of the filter width 7 and fixed uo = 0. 
The blue, green, red, and yellow areas refer to values of 
7T0 = 2000, 40, 10, and 2, respectively. The other control 
parameters are fixed as R = 0.7 and r = To/2. The parame- 
ters of the system are as in Fig. [1] 



shifted up towards larger values of R. The left boundary, 
which corresponds to if m i n , is shifted towards larger val- 
ues of K. This effect is very small for 7T0 = 40 and 10. 
For 7T0 = 2, the effect is much more pronounced. 
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FIG. 15: (Color online) Domain of control in the (r, K) plane 
for different values of the filter width 7 and fixed uoo = 0. 
Panels (a), (b), (c), and (d) correspond to values of 7T0 = 
2000, 40, 10, and 2, respectively, where the time delay is fixed 
at r = To/2 and R = 0.7. The parameters of the system are 
as in Fig. [1] 

Figure [15] shows the domain of control in the (if, r) 
plane for different values of the filter width 7 and fixed 
uoq = 0. Panels (a), (b), (c), and (d) correspond to values 
of 7T0 = 2000, 40, 10, and 2, respectively, where the 
time delay is fixed at r = To/2 and R = 0.7. The case of 
7T0 = 2000 looks very similar to the case of the unfiltered 
system, as shown in Fig. [2] For decreasing value of 7, 
the tongues shrink in K direction. The minimum value 
of K increases and the maximum value becomes smaller. 
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Additionally, the tongues are bent down towards smaller 
values of the time delay r with decreasing 7. This leads 
to an optimal r that is smaller than in the unflltered case. 




2 4 6 2 4 6 



KT KT 

FIG. 16: (Color online) Domain of control in the (r, K) plane 
for different values of the filter's mean value uo and fixed 
7T0 = 10. Panels (a), (b), (c), and (d) correspond to values 
of ujqTq — 0, 7r, 27T, and 47T, respectively, where the time 
delay is fixed at r = To/2 and R = 0.7. The parameters of 
the system are as in Fig. [T] 

In Fig. EH the filter width is fixed to 7T0 = 10, the do- 
main of control is shown for different values of the filter's 
mean values ujq. Panels (a), (b), (c), and (d) correspond 
to values of ujqTq = 0, 7r, 27r, and 47r, respectively, where 
the time delay is fixed at r = To/2 and R = 0.7. The size 
of the domain of control is only slightly changed with in- 
creasing uoq. The domains are flattened on the upper side 
for larger values of uoq. The region of optimal control, de- 
noted by yellow color code, is additionally shifted slightly 



towards larger values of the feedback gain K and shrinks 
in r direction. Overall, the variation of uoq has very little 
effect on the domain of control in the (if, r)-plane. 



VI. CONCLUSION 

In conclusion, we have shown that extended time- 
delayed feedback can be used to stabilize unstable steady 
states of focus type. By introduction of an additional 
memory parameter, this method is able to control a larger 
range of unstable fixed points compared to the original 
TDAS scheme. However, the degree of stability, mea- 
sured by the absolute value of the real part of the eigen- 
value, is generally decreased. We have investigated the 
domain of control in various one- and two-dimensional 
projections of the space spanned by the control param- 
eters. Furthermore, we have discussed also effects of 
nondiagonal coupling, nonzero control loop latency, and 
bandpass filtering of the control signal, which are rele- 
vant in an experimental realization of the ETDAS con- 
trol method. We found that a proper adjustment of the 
time delay is able to compensate for reducing stabiliza- 
tion abilities of the control method, for instance, due to 
latency or feedback phase. 

We point out that the results obtained in this paper 
are accessible for applications in the context of all-optical 
control of intensity oscillations of semiconductor lasers as 
investigated in [37[. 
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